The Spectral Energy Distribution of Galaxies 
Proceedings lAU Symposium No. 284, 2011 
R.J. Tuffs & C.C.Popescu, eds. 



©2011 International Astronomical Union 
DOI: 00.0000/XOOOOOOOOOOOOOOOX 



Constraints on the Optical-IR extragalactic 
background from 7-ray absorption studies 

Luigi Costamante^'^ 

^ Stanford University, USA & INAF-OAB, Milano, Italy 
Now at: Dept of Physics, University of Perugia, Perugia, Italy 
email: luigic2011(9gniail . com 

■ Abstract. Very high energy (VHE, > 0.1 TeV) gamma-rays from extragalactic sources, inter- 
acting by 7 — 7 collisions with diffuse intergalactic radiation fields, provide an alternative way 
to constrain the diffuse background light, completely independent of direct measurements. The 
, ^ limits depend however on our knowledge of the physics of the gamma-ray sources. After clar- 
O ^ ifying the interplay between background light and VHE spectra, I summarize the extent and 
validity of the obtainable limits, and where future improvements can be expected. 
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1. Introduction 

VHE 7-rayf[3 from extragalactic sources provide an alternative and completely indepen- 
rS i' dent way, with respect to direct measurements, to probe the diffuse Extragalactic Back- 
ground Light (EBL, from UV to far infrared wavelengths; see e.g. lHauser fc Dwek 200l|) . 
This approach is based on the study of "absorption" features imprinted on the GeV-TeV 
^ spectra by the interaction of 7-rays with EBL photons through the pair-creation process 
■ (77 e+e"; see e.g. lAharonian 20011 and refs therein). Blazars represent a very useful 
class of 7-ray beamers, being numerous over a wide range of redshifts, very luminous 
(enhanced by relativistic beaming) and long-lasting sources. However, they are far from 
I> . being standard candles, and therefore the constraints on the EBL are always heavily 
' dependent on our understanding of the blazar intrinsic emission and physical properties. 
. Nevertheless, it is possible to derive meaningful constraints if the intrinsic blazar prop- 
erties implied by certain EBL spectra are very far from the known or expected range, 
and/or if they become inconsistent with other aspects of the blazars spectral energy 
distribution (SED). Observations in the last 5 years allowed a substantial progress, pro- 
. viding the stringest constraints to date. Here both the main results and the limits of 
their validity are discussed. 



2. Diagnostics: how absorption deforms TeV spectra 



^ _ The convolution of the 7-7 cross section with the EBL spectrum yields an energy- 
" " " dependent attenuation. Because the resulting optical depth r mostly increases with 7-ray 
energy, the observed 7-ray spectrum emerging at the end of the travel path is steeper 
than the initial, intrinsic spectrum emitted by the source (namely, the photon index 

f In the following. 'VHE' and 'TeV will be used as interchangeable terms, meaning about 
±1 decade around 1 TeV. The same for 'HE' (high-energy, >0.1 GeV) and 'GeV terms. In fact, 
TeV and GeV are the energies at which atmospheric Cherenkov telescopes and the Fermi-LAT 
detector have the best sensitivity and resolution, respectively, for most astronomical spectra. 
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Figure 1. Left: SED of the EBL, as was discussed in year 2005 (fromlAh aronian ct al. 2006|) . 
Open symbols show lower limits from galaxy counts, while filled symbols correspond to direct 
estimates. On the upper axis it is plotted the TeV energy corresponding to the peak of the 7 — 7 
cross section. The two curves (identical except for normalization) are drawn as simple shapes 
reproducing the main EBL features. Right: the attenuation curves corresponding to the 
EBL curves in the left panel (lower solid and dotted lines), and at two different redshifts for 
the upper EBL curve (dark and grey solid lines). The attenuation curves represent directly the 
shape of an observed VHE spectrum, if the initial spectrum were a power-law parallel to the 
upper axis. All calculations performed assuming Ho = 70 km/s/Mpc, ilm = 0.3, = 0.7. 



^obs ^ Tint). This can be seen in Fig. [1] where the attenuation factor e~'^ is plotted as 
a function of energy. Limits on the EBL can thus be derived if the intrinsic spectrum 
required by a particular EBL intensity (or shape) is too anomalous -typically too hard- 
with respect to the knovifn physics of the source (e.g. ICostamante et al. 2004|) . Several 
points can be noted from Fig. [1] 

1) different EBL wavebands (having different slopes) affect each part of the 7-ray 
spectrum in a different way, yielding a complex deformation shape. However, over some 
ranges (namely 0.2-1 TeV and 2-8 TeV), it approximates a power-law shape. That is, 
if the intrinsic spectrum is a power-law, the observed spectrum can be well fitted by a 
power-law of steeper index. 

2) The amount of steepening AL increases both with distance and EBL intensity. 
Thus redshift gives leverage: the same uncertainty in EBL flux AFebl causes a higher 
AL at larger distances. More distant sources provides therefore more sensitivity for EBL 
constraints, but at the same time they suffer stronger attenuation, yielding much lower 
statistics. As a result, there is an optimal range of redshifts that provides the best com- 
promise between EBL-sensitivity and statistical error, and thus gives the most stringent 
EBL limits. With the present detectors, it is around z 0.15 — 0.25. 

3) If the EBL number density n{e) oc e^*^ (i.e. the EBL SED oc A^^^)^ ^^le optical 
depth t{Ej) becom es cx E^~^. Whe re /3 « 1, the optical depth becomes therefore inde- 
pendent of energy (|Aharonian 200l|) . In such case there is no steepening: the absorbed 
spectrum reproduces the original shape, simply attenuated by a constant factor. This is 
what partly happens in the EBL waveband 3-10 /xm, according to all most recent EBL 
calculations (e.g. IFranceschini et al. 20081 |Dominguez et al. 2011| IGilmore et al. 20111 
IKneiske fc Dole 2010p . It causes a flattening feature in the attenuation curve between 
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^1-2 and 8-9 TeV, where the observed 7-ray spectrum partly recovers its intrinsic slope 
and thus can appear harder than in the 0.2-1 TeV bandQ 

4) Therefore, NO cutoff is produced by EBL absorption between 0.2 and 8 TeV \ Any 
cutoff or steepening seen in this energy range should be intrinsic to the source spectrum. 
Fig. [T] also clarifies why almost all "data points" reported on the so-called "gamma-ray 
horizon" plots in literature (showing the at which r = 1 as a function of redshift) are 
basically meaningless: there is no way to measure an EBL cutoff between 0.1 and 10 TeV 
using only VHE data. The only two EBL cutoffs are between the HE and VHE bands, 
i.e. between 20 and 200 GeV (thus requiring a Fermi-LAT spectrum), or above 8-10 TeV 
(for which only two blazars have been detected with sufficient statistics so far). 

5) The 7-ray steepening is determined by the difference in optical depth between 
the two ends of an observed VHE band. Therefore any spectral steepening (and conse- 
quently any EBL upper limit) can always be counteracted and canceled by an appro- 
priate increase of the optical depth (i.e. of the EBL flux) at the lower 7-ray energies, 
so to equalize the optical depths. In the 0.1-1 TeV range, this is obtained by increasing 
the UV flux with respect to the NIR flux (see lAharonian et al. 2006|) . Therefore, any 
limit on the EBL always depends on the assumption on the UV flux, by definition. All 
claims of "new" methods for "model-independent" or "shape-independent" EBL limits 
(e.g. IMazin fc Raue 20071 IHinton fc Hofmann 2009[) are therefore essentially incorrect, 
even if their (often hidden) assumptions are reasonable and in fact commonly used (e.g. 
assuming that the UV flux is always lower than the Opt-NIR flux). 



3. Optical-NIR constraints 

Until 2005, the large uncertainty especially around 0.8-3 /im (a factor 5-10 x between 
lower limits and direct estimates) caused a fundamental ambiguity in the interpretation 
of blazars gamma-ray spectra. An observed TeV spectrum could be the result either of a 
hard intrinsic spectrum attenuated by a high-density EBL, or of a soft intrinsic spectrum 
less absorbed by a low-density EBL. In both cases, the required intrinsic TeV spectrum 
was well within the typical range of spectra shown so far by blazars, so no useful constraint 
could be derived (with the exception of the 2-10 /im range, see lAharonian et al. 2003| . 
This ambiguity essentially undermined our capability to study blazars, since it corre- 
sponded to large differences both in energy and luminosity of their SED 7-ray peak. 

In 2005, however, the HESS observations of the BL Lac objects lES 1101-232 (z=0.186) 
and H 2356-309 (z=0.165) provided a fundamental breakthrough (jAharonian et al. 2006|) . 
The observed 7-ray spectra (detected between 0.2 and 1-3 TeV) were much harder than 
expected for their redshift, implying extremely hard intrinsic spectra (Tint ^ 0) in case 
of high EBL densities (Fig. [2|). Such hard spectra were never seen in the closer, less ab- 
sorbed objects and were at odds with all the known blazar physics and phenomenology. 
They were also not supported by same-epoch multiwavelength observations of their SED 
(jAharonian et al.~2007a|) . A low EBL intensity could instead accomodate all the new re- 
sults within the typical range of blazar properties. Assuming that the intrinsic spectra 
were not harder than r=1.5 an upper limit could be derived, which resulted very close to 
the lower limits given by the integrated light of galaxies (details in lAharonian et al. 2006p . 
This result is now further corroborated by several other sources and observations (e.g. 
lAharonian et al. 2007b1 lAcciari et al. 20091 [Albert et al. 2009|) . 

f In fact, 7-7 absorption can make a gamma-ray spectrum even /larder than originally emitted, 
when /3 < 1 (e.g. fo r a planckian distribut ion). This is indeed an effective way to explain very 
hard 7-ray spectra (|Aharonian et al. 20081) . 
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Figure 2. Left: The HESS spectra of lES 1101-232, corrected for absorption with three dif- 
ferent EBL SEDs, as labelled f from [Xharonian et al. 2006|) . Lower points (red): observed data. 
Upper points (blue): absorption-corrected data. The lines show the best fit power-laws to the 
reconstructed spectrum, and the corresponding shape after absorption. Right: upper limit for 
the EBL (black marked region), from the assumption of a blazar spectrum F ^ 1.5. 

This result has three niain consequences: 1) it pins down the origin of the EBL, show- 
ing that at these wavelenghts it is strongly dominated by the direct starlight from 
galaxies, excluding a strong contribution from other sources like Pop-Ill stars (e.g. 
ISantos et al. 2002)) : 2) it means that the intergalactic space is more transparent to 7- 
rays than previously thought, thus enlarging the 7-ray horizon; 3) it strongly reduces the 
ambiguity on blazars TeV spectra. 

3.1. The (in)famous limit of T — 1.5. 

There is a lot of confusion in literature on the validity of the r=1.5 hardness limit 
and on the reasons of its adoption. It is important to understand that this is NOT the 
hardest possible theoretical spectrum in blazars, and was never introduced as such (see 
lAharonian et al. 2006p . Indeed there are many possible mechanisms to produce extremely 
hard gamma-ray spectra in blazars, such as bulk- motion comptonization ([Aharonian 200 1|) . 
internal absoption on narrow-banded photon fields ([Aharonian et al. 2008p . uncooled 
particle acceleration spectra or fine-tuned shock acceleration (yielding however T ~ 1.2, 
e.g. IStecker et al. 2007)) . a low-energy cutoff in the particle distribution at very high 
energy (e.g. jKatarzynski et al. 2006 ILefa et al. 2011)) . or relativistic Maxwellian distri- 
butions (e.g. jSauge' fc Henri 2004[ILefa et al. 2011|) . 

The r=1.5 limit is a reference value, a benchmark that so far marks the borderline be- 
tween reality and speculation. Observations do show that blazars can have photon spectra 
with r > 1.5 (directly seen in synchrotron or inverse Compton emissions), and these can 
be produced with standard shock acceleration and cooling mechanisms without invoking 
special conditions or fine-tuning. On the other hand, much harder spectra have never 
been significantly detected in blazars so far, at high particle energies, neither directly or 
by synchrotron emission. The alternative scenarios, though not excluded by observations, 
are not suggested by them either, and in some cases require extreme fine-tuning to avoid 
confiicts with multiwavelength SED data. 

Furthermore, the r=1.5 value should not be considered a sharp, hard limit. At any 
given redshift, EBL absorption establishes a one-to-one relation between the observed 
and intrinsic 7-ray spectral slope. Thus a different assumption shifts the EBL limit 
accordingly, but only a large change (e.g. Tint « 1) can effectively impact the conclusion 
of a low EBL. For example, in the case of lES 1101-232, a shift of AF = 0.3 (e.g. by 
assuming Tint ^ 1-2) affects the EBL limit by no more than 10% (|Ah aronian et al. 2006)) . 
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absorbed r,„,= 1.5 (fit 0.2-1 TeV) 




Figure 3. Left: the SED of lES 0229+200, showing the correct X-ray data from Swift (fiUed 
grccn/grey points), at much lower flux than the wrong X-ray data (open triangles) used in 
[Tavecchio et al. 20091 Right: observed VHE photon indices of all detected TeV blazars as a 
function of redshift. The lines represent the expected observed (i.e. absorbed) photon index of 
a source with an intrinsic spectrum F = 1.5, when fitted in the range 0.2-1 TeV, for different 
EBL levels (as labelled). The separation of the data from the lines at each redshift shows how 
much the intrinsic spectrum is harder or softer than 1.5. Adapted from Costamante 2006. 

The Opt-NIR upper limit is therefore robust to such changes in blazar assumptions, which 
are of the same order of the statistical and systematic uncertainty in the 7-ray spectra. 

The only relevant issue left, therefore, is to assess -theoretically and observationally- 
if blazars do have emission components characterized by extremely hard spectra (e.g. 
r = 0.7 or less) or not, and when they produce them. 



4. Why a low EBL density is still the most likely solution. 

Despite the several possible scenarios for producing very hard spectra in blazars, there 
are two main reasons why a low EBL density seems still the correct solution. 

The first concerns the emitting particles spectrum. Most of these scenarios invoke very 
hard particle spectra, either as a pile-up/Maxwellian distribution or as a low-energy cut- 
off, and suppressing radiative cooling. These mechanisms look feasable even in a one-zone 
synchrotron self-Compton scenario for blazars (see lLefa et al. 20111 for a recent compre- 
hensive discussion). However, since the synchrotron emission (which in these objects 
correspond to optical-X-ray energies) traces directly the particle spectrum, such hard 
features should appear in the synchrotron spectrum as well, at some energies and at 
least in some epochs. Instead, they have never been observed in more than 30 years of 
X-rays and optical/ir observations of blazars, requiring a " cosmic conspiracy" to always 
hide the hard emission below a more normal component or in seldom-observed bands 
(mm ?). Intriguingly, there was a recent claim of possible observational evidence for a 
low-energy cutoff in the Swift data of lES 0229+200, as an unusually high X-ray-to-UV 
flux ratio ([Tavecchio et al. 2009p . However, a bug was found in the calculation of the 
effective area of the X-ray data, and the correct result (see Fig. [3]) show a spectrum 
in line with all other objects and observations. Future observations with ALMA (mm) 
and NuSTAR (hard X-ray) may provide further insights. Dropping the requirement of a 
hard particle spectrum, internal absorption on a narrow-banded photon field can make 
the gamma-ray spectrum extremely hard (jAharonian et al. 2008^ However, large fiuxes 
in the GeV range are required, and these seem now excluded by the low fluxes or non- 
detection of the hard- TeV sources with Fermi-LAT (jAbdo et al. 2010p . 



6 



Luigi Costamante 



MKN 501 z=0.034 HEGRA Dominguez 201 1 EBL 




X[/im] Log E [TeV] 



Figure 4. Left: SED of the EBL with the most recent models. The grey band show the EBL 
upper limits in the Opt-NIR range (with the +la statistical uncertainty from the VHE data) 
re-calculated using the EBL shape by Franceschini et al. (2008). Right: HEGRA spectrum of 
Mkn 501, as observed (red, lower points) and absorption-corrected with two recent EBL-model 
calculations ( [Dominguez et al. 2011| and IGilmore et al. 20111 upper black-blue and grey-green 
points, respectively). Both cases show an upturn. 

The second reason for a low EBL density (and thus for Tint ^ 1-5) is that a high 
EBL would require a dramatic change of blazar properties in a very narrow range of 
redshifts (see Fig. [3]). A high EBL would create a sharp dicothomy in blazars spectra 
around z = 0.15: all sources below this redshift would have F ^ 1.5, while all blazars 
above this redshift would have always F < 1.5, all other properties being equal or very 
similar. There is no known reason, observational bias or evolution parameter that can 
explain such an abrupt change of properties in such small redshift range. Instead, a low 
EBL intensity makes the range of TeV blazar spectra consistent among each other at any 
redshift sampled so far. 

Therefore, though not yet "bullet-proof, until proven otherwise a low EBL seems 
the preferable solution according to Occam's razor, (see also lMadau fc Silk 20051 for the 
cosmological problems to explain a high EBL flux) . 

5. The Chain of Constraints 

Since the EBL photon field is the same for all sources, and evolving in redshift, the 
spectral-hardness limits from different blazars and in different energy bands can be com- 
bined to constrain the EBL over a wide range of wavelengths, and in a more stringent way 
than allowed by each single object ([Costamante et al. 20041 IDwek fc Krennrich 20051 
IMazin fc Raue 2007p . Each new spectrum can take advantage of the previous limits, 
forming a chain of constraints that starts with (and depends on) the assumption on the 
UV flux. With the reasonable hypothesis that the UV background is lower than the Opt- 
NIR one (e.g. around 2-6 nW/m2 sr), blazars spectra between 0.1 and 1 TeV pin down 
the Opt-NIR flux very close to the galaxy counts limit ([Aharonian et al. 2006^ . 

With this limit, the blazar spectra between - 1 and 10 TeV (lES 0229-H200 [Aharonian et al. 2007cl 
and lES 1426-I-428. Aharonian et al. 2003^ constrain the EBL spectrum between 2 and 
10 /j,m to a slope close to A~^, again very close to the lower limits from galaxy counts 
(Aharoni an et al. 2007cp . This is possible because the Opt-NIR limit precludes the pos- 
sibility of increasing the flux at '^1 /Ltm to reduce the difference in optical depth between 
1 and 10 TeV (see Sect. 2). 



EBL constraints from gamma-ray absorption studies 



7 



The consequent upper limit at 10 /xm constrains the rising of the EBL SED towards 
the far-infrared hump, which is determined by warm and cold dust emission. This band 
is probed by 7-rays between ~ 8 and 100 TeV. So far, only two close-by objects have 
been measured up to 15-20 TeV, thanks to HEGRA observations: Mkn 421 (2:=0.031, 
lAharonian et al. 1999)) and Mkn 501 ^2:=0.034. rAharonian et al. 200l|) . Interestingly, the 
HEGRA spectrum of Mkn 501 seems to have problems with the most recent EBL calcu- 
lations, which cause an up-turn or pile-up at the highest energies (see Fig.|4]). This is the 
modern version of the so called "TeV-FIR background crisis" (see e.g. lAharonian 200ip . 
caused years ago by a first very high estimate of the 50 /im EBL flux. However the in- 
formation on warm dust and the 7-ray statistic are still insufficient to draw any sensible 
conclusion. More data > 10 TeV are needed (hopefully from CTA). 



6. Future and Conclusions 

VHE data sample mainly the local EBL (z = — 0.5). To probe the EBL evolution 
over cosmic time and in the UV band, sources up to 2; = 3 — 4 and data in the 10- 
100 GeV range are needed. The first results from Fermi-LAT (jAbdo et al. 2010bl) are 
in agreement with the VHE limits and with the most recent EBL evolution models. 
To this respect, however, GRBs might turn out to be more useful 7-ray sources, given 
that EBL and blazar-intrinsic evolutionary effects might be very difficult to disentangle 
(jReimer A. 2007p . 

At VHE, instead, CTA is NOT expected to improve the Opt-NIR limits per se, contrary 
to various claims. The reason is twofold. First, there is not much room for improvement 
! The present upper limits already match the lower limits by galaxy counts quite closely 
(see Fig.S]) Second, even with infinite 7-ray statistics, it remains the unavoidable system- 
atic uncertainty in blazar modeling. The small AF ~ 0.1 — 0.3 induced by the residual 
EBL uncertainty between lower and upper limits can be typically accomodated with very 
small changes in blazar parameters, for a given set of multiwavelength data. 

The main question to be addressed is if our blazar assumptions are correct (and in 
such case the EBL is already pinned down) or are not (and in this case by how much) . In 
the latter case, the entire construction must be revisited, since it is our understanding of 
the gamma-ray sources that changes completely. It is on this aspect, on blazar physics, 
that CTA is expected to provide the most significant improvements and to test our 
assumptions by finding counter-examples. This can be done by monitoring more low- 
rcdshift objects in the TeV range (to detect directly F < 1) and by measuring more 
high-redshift spectra (to find those for which F — 1.5 is incompatible with galaxy counts; 
ICostamante 2007p . 

In conclusion, a low EBL close to galaxy counts seems the "convergent solution", 
despite some uncertainties on blazar physics. However, there are still some fundamental 
aspects of blazars' acceleration and emission mechanisms not yet understood, which 
present and future observations are testing. The EBL limits obtained so far are robust 
against small changes in blazars assumptions, but they all depend, by construction, on 
the assumption of a low diffuse UV flux. 
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DISCUSSION 

Q (K. Mattila): Does MAGIC give additional (or better) constraints to the EBL at 
optical wavelengths, where it should be more sensitive than HESS because of its lower 
7-ray energy range ? (I am referring to the 2008 Science paper on 3C279). 

A: In this case (3C 279), not really. The problem is the much higher uncertainty in the 
spectrum. This is also caused by 1) the smaller detected band (less than half a decade, 
with respect to a full decade in energy for lES 1101-232 and the other BL Lacs), and 2) the 
much lower S/N. The lower 7-ray energy range does probe slightly shorter wavelengths, 
but it also reduces the sensitivity of the 7-ray spectrum to EBL changes, with respect to 
the 0.2-2 TeV range, despite the higher redshift. Note that in the blue-UV range ("max 
EBL" curve in their Fig. S2), the given values are basically arbitrary, and should not 
be considered limits at all (a higher UV flux makes also the Magic spectrum softer, see 
Sect 2). In summary, that observation does not provide more stringent constraints on 
the optical background, though the detection itself corroborates the idea of a low overall 
EBL (adopting historical GeV peak fluxes for 3C279). It is also important to remind 
that, for this source, there is the further systematic uncertainty of internal absorption 
on UV photons from the broad line region (this object has broad lines in the optical 
spectrum) . 



